Localizing gravity on Maxwell gauged CP 1 model in six dimensions 
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We shall consider about a 3-brane embedded in six-dimensional space-time with a negative 
bulk cosmological constant. The 3-brane is constructed by a topological soliton solution living 
in two-dimensional axially symmetric transverse subspace. Similar to most previous works of six- 
dimensional soliton models, our Maxwell gauged CP 1 brane model can also achieve to localize 
gravity around the 3-brane. The CP 1 field is described by a scalar doublet and derived from 0(3) 
sigma model by projecting it onto two-dimensional complex space. In that sense, our framework is 
more effective than other solitonic brane models concerning with gauge theory. We shall also discuss 
about linear stability analysis for our new model by fluctuating all fields. 

PACS numbers: ll.10.Kk, 11.27.+d, 04.50.-h 



I. INTRODUCTION 

In recent years, manifolds living in a higher dimen- 
sional space-time well known as D-branes have fascinated 
numerous physicists. They have been derived from topo- 
logical soliton solutions in string theory [![, which arise 
naturally in ten-dimensional supergravity or string/M- 
theory 0, 11,1,0. It is an old idea that space-time may 
have more than four dimensions and the extra ones are 
unobservable for ordinary energy scales @, H, H OH HH • 
There are generally three different possibilities of extra- 
dimensional scenario: extra-dimensional space is com- 
pact [ljj , non-compact but has a finite size [H [3, EH > 
and non-compact with an infinite size In particular, 
Ref. [l(| leads to localized gravity around a 3-brane with 
non-zero tension in anti-de Sitter (AdS) space though the 
extra-dimensional space with an infinite size. 

The brane proposed in Refs.[l5|, [l6[ and its general- 
izations, i.e., models of Randall- Sundrum (RS) type, are 
essentially static point-like external sources in the extra 
dimensions. On the other hand, an increasing interest re- 
cently has focused on study of gravitating thick (or fat) 
defects embedded in higher dimensional space-time with 
codimension one or more [13, H O H Hi El El ■ Our 
main proposal is to construct a new 3-brane which is 
described by special classical solutions, i.e., topological 
solitons of a field-theoretical Lagrangian. This approach 
is inspired by the D-branes where the solutions are con- 
structed by topological solitons in string theory. In that 
sense, the solitonic branes may be more natural rather 
than another descriptions by the delta-function like po- 
tential. 

3-branes in relevant previous works are constructed by 
the classical solutions in gauge theory; such as kink in 
five dimensions (5D), Abelian-Higgs vortex in 6D, and 
t'Hooft-Polyakov monopole in 7D 0, El El, El HI- 
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In particular, codimcnsion-2 brancworld models con- 
structed from Abelian-Higgs vortex in 6D are studied for 
a only case of flat 3-brane but also curved 3-brane [27| . 
On the other hand, our 3-brane is written by a Maxwell 
gauged CP 1 model whose origin is O (3) sig ma model in 
(2+l)-dimcnsional field theory [H M, El El El El . 
0(3) sigma model appears in various aspects of physics, 
and has richer topological classes than Abelian-Higgs vor- 
tex; lump, gauged vortex, baby skyrmion, and instanton. 
The 0(3) sigma model can possess finite energy soliton 
solutions by adding some additional terms or gauges into 
the Lagrangian. The solitons in gauge theory are strongly 
restricted by the configurations of their gauge groups; on 
the contrary, our framework is flexible and thus can be 
easily applied for various background geometry 

Another aim of the present paper is to analyze lin- 
ear stability of our new solutions by fluctuating all fields 
[35L EH, [!?], EH . Study of a linear stability and the sec- 
ond order correction of gravity and coupled fields for 
models of RS t ype have been done by numerous au- 
thors [H SO, 5111 51 51, 51, 51- Also, analysis for 
gravitating thick defects embedded in higher dimensions 
are found in the literature; for 5D [H 53, 5H, and for 
6D [H, H3, |H HI. (Note that the model in Ref. HI 
is constructed by gravitating multidefects in five dimen- 
sions). The studies for thick defects, however, are works 
in progress since the topological defects used in the lit- 
erature are quite complicated structures. More seriously, 
stabilized mechanisms for the models of RS type, like 
"moduli stabilization" , have not been found. Analyz- 
ing the linear stability of our obtained solutions is thus 
worthwhile to tackle. 

The plan of this paper is as follows. In Scc[TT]wc build 
up the model and the resulting equations of motion. In 
Scc lIIII we describe the boundary conditions for the mat- 
ter fields and the warp factors. In Sec lIVI we mainly dis- 
cuss about the asymptotic behavior of the solutions of the 
model. In Sec|V]we introduce the methods for solving our 
boundary-value problem. Some typical results localizing 
gravity around the 3-brane are shown in this section. In 
Scc lVII wc give a detailed analysis for the stability of our 



brancs for linearized gravitational perturbations. Finally, 
in Sec I VIII some conclusions and summary of this paper 
are drawn. 



II. SIX-DIMENSIONAL MODEL AND FIELD 
EQUATIONS 

In this section we shall construct Maxwell gauged CP 1 
model combining with general relativity in six dimen- 
sions. As in the previous works for braneworld scenarios, 
an action of the models is written as a coupled system 
with gravity and solitons. The total action of our six- 
dimensional model is written as 



S = Shi 



where Sbrane ^ s Maxwell gauged CP 1 model action and 
5 g rav is the six-dimensional generalization of Einstein- 
Hilbert gravity. The explicit form of Einstein-Hilbcrt 
gravity is given by 



A, 



(2.1) 



where R is the six-dimensional Ricci scalar, is the 
bulk cosmological constant, \ = 8ttGq = %-k/Mq and M 6 
denotes the six-dimensional Planck mass. On the other 
hand, the brane action 1 is of the form 
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k(V M z a yV M z a -fi 1 - z a * (a 3 ) ab z b 



(2.2) 



where a 3 is the third component of the Pauli matrix. 
The action (|2.2[) is strongly motivated by the (2+1)- 
dimensional field-theoretical model [U [33j|. The cou- 
pling constants g, k, [i describe those of the field strength 
of the gauge field, the kinetic term of CP 1 and the po- 
tential term. Their dimensions arc [M~ 2 ], [M 4 ], [M 6 ] in 
the n.u.(M denotes a unit of mass), respectively. The 
gauge covariant derivative T>m in Eq. (|2.2p is defined in 
terms of the (7(1) gauge field Am such as 



connection substantially works as a U(l) gauge field [30| . 
Thus we replace the composite field connection with the 
U(l) gauge connection Am and treat as the independent 
field variable [33j|. It should be noted that our model 
with this replacement is not to be completely equivalent 
to 0(3) sigma model. 

Formally by varying the actions (|2.ip and (|2.2p with 
respect to the field variables, one can obtain the classical 
equations of motion 
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(2.3) 



Rain — -^GainR — x (Tmn — ^Gmn) ■ 
V M {V M z a ) + iA M V M z a + 1<jfz b = 0, (2.4) 



|JV I 1 2 . ... ( ac>N . 



V M P = 8kgA N \z\ 2 + iikg (z a d N z a * - z a *d N z a ) , 

(2.5) 

where Vm is the covariant derivative with respect to the 
metric tensor and 
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is the six-dimensional energy-momentum tensor. 

At the present paper, we would like to consider about 
a warped six-dimensional space-time with axially sym- 
metric two extra dimensions. The ansatz is imposed on 
the metric tensor 



ds z = Gmn<1x m dx N 
= M 2 (r)rj )iV dx fl dx' / 



+ dr 2 + L 2 (r)d6 2 1 (2.7) 



where r and 9 are, respectively, the bulk radius and the 
bulk angle, ?7 M „ denotes four-dimensional Minkowski met- 
ric tensor and the function M(r),L(r) are often called 
warp factors in brane world scenarios. For the matter 
fields, we explore the solutions for the CP 1 doublet field 
and the U(l) gauge field of the form: 



z - 
A„ 



{zi,z 2 ) T 
-0, A r 
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(2.8) 
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In standard CP 1 model, the Am is not an independent 
field but is the composite field connection in terms of CP 1 
field defined as Am = iz a *dMZ a - The composite field 



1 The conventions for indices of the present paper are the follow- 
ing: the capital Latin indices run about from to 5, the Greek 
indices run about from to 3, and the small Latin indices which 
denote the complex projection (CP) field, therefore the small 
Latin indices run about from 1 to 2. 



where n is the winding number of the gauge field and the 
CP 1 doublet z satisfies the constraint z^ z = 1. 

The classical equations of motion of the present system 
arc then 



f" + (4m + l) f - 2 (J§ + 7 sin f) = 0, 

5" + (4m - i)a! - v a = 0, 

£' + 3m' + I 2 + 6m 2 + 3m^ = a (r - (5) 

4m + 10m 2 = a (rg - /3) , 



4m^ + 6m 



a tv 



(2.9) 

(2.10) 
(2.11) 
(2.12) 
(2.13) 
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where a := \k is the dimensionlcss gravitational coupling 
constant, (3 := Ai,/k 2 g is the dimensionlcss bulk cosmo- 
logical constant, and 7 := fi/k 2 g is the dimensionlcss 
coupling constant. Besides the dimensionless coordinate 



kgr 



is introduced and the prime in the equations denotes the 
derivative with respect to this dimensionless coordinate. 
The matter fields f(r), a(r) and the warp factors M(r), 
L(r) are rewritten in terms of the new coordinate, that 
is, 



a(r), 

= v^ffi(r). 



f(p) ■= f(r), HP) ■■ 

M(p):=M(r), C( P ) 



The functions m(p) and £(p) arc defined by the above 
functions as 



dlnM(p) , . d\nC(p) 
™{P) := , £{p) ~ 



dp 



dp 



(2.14) 



which seem to have some similarity with the Colc-Hopf 
transformation appearing in the integrable theory. The 
components of the energy-momentum tensor Tj^ go to 
be dimensionlcss ones 

•»«=-? —if- p-Jr-id-cos/), 
7*00 : =T = -4 /,2 + 2^ + ^- 7(1 - cos/) ' 

(2.15) 



III. BOUNDARY CONDITIONS 
A. Soliton 

Existence of topological soliton solutions is inferred 
from Derrick's scaling argument [53| in which, if soliton 
exists, a stationary point of the energy in the field config- 
uration should be stationary against all variations includ- 
ing spatial rescaling. Also, the soliton solutions always 
have the lower energy bound which is defined by their 
topology. For the U(l) gauged CP 1 model with Maxwell 
or Chern-Simon term, such existence proof was confirmed 
numerically [3l|, [32|, [33j]. Besides the original Maxwell 
gauged CP 1 model in (2+1) dimensions can have topo- 
logical soliton solutions [331 ] ■ These analysis clearly in- 
dicate that the solitons can be stabilized without any 
higher order terms, like the Skyrme term in (3+1) di- 
mensions [29|. 

If we apply these models to the six-dimensional space- 
time, both the CP 1 field and the U(l) gauge field should 
go to zero at infinity and be regular at the origin in the 
extra-dimensional space, which exactly agree with the 
topological requirements of the model. These conditions 
are 



/(0) = 7T, 

5(0) = n, 



lim f(p) = 0, 

p — >OC 

lim a(p) = 0. 

p — >oo 



(3.1) 



As is well known that <CP N model has a close relation to 
0(3) sigma model. For the case oi N = 1, one can easily 
obtain the 0(3) sigma model by using the transformation 



where the dimensionlcss quantities v,v a ,Vf have been de- 
fined as 



v := a 2 + n(n — 2a) sin 2 — , 



dv 
da 



2 / 

2 \ a — n sin 



2 / ' 



dv n(n — 2a) 
Vf : 77/ 2 Sin/ - 



(2.16) 



Before beginning our analysis, let us consider basic 
properties of the field equations (|2.9p - (|2.13p . Clearly 
Eqs. (|2.9p and (|2.10|) are the dynamical equations because 
they are second order differential equations of the fields 
f(p),a(p). Since Eqs.( |2~TTj) and (|2~T2|) contain £', to', i.e., 
the second derivative of the warp factors, these equations 
are the dynamical equations, too. Eq. (|2.13p is comprised 
of the first derivatives only. Thus it means that the equa- 
tion works constraint equation for the dynamical 
fields f(p),a(p),A4(p),C(p). As a result, we must treat 
a numerical problem of a series of four dynamical equa- 
tions with one constraint equation. 



where 4> a , a — 1,2,3 are scalar triplet fields, z a are dou- 
blet fields, and <r a are the three Pauli matrices. In the 
sense, the original CP 1 model is essentially 0(3) sigma 
model when it is written in terms of the 4> a . Moreover, in 
the view of 0(3) sigma model, the topological boundary 
condition is interpreted as the south pole configuration 
at the origin and the north pole configuration at infin- 
ity. The details of above discussion can be seen, e.g., in 
Ref . [30l] . As mentioned before, our model however is not 
equivalent to 0(3) sigma model completely. 



B. Geometry 

Regular geometry at the origin defines the bound- 
ary conditions for the warp factors. An effective ac- 
tion derived from six-dimensional gravity action gener- 
ally has two kinds of singularity at the origin, which are 
called conical and curvature singularity problem (e.g., in 
Ref. [54| ) . The boundary conditions for geometry are in- 
troduced in order to exclude these serious difficulties in 
the six-dimensional model. In this paper they arc given 



3 



by 



M'(0) = 0, C(0) = 0, £'(0) = 1. 



(3.2) 



One simply fixes M(Q) = 1 since the value of M(0) is an 
arbitrary constant. The boundary conditions for £ and 
M' correspond to the conical and the curvature singu- 
larity problem at the origin, respectively. Also, we con- 
centrate on the problem of the regular geometry without 
the conical singularity so that we employ the boundary 
condition £'(0) = 1. 

In order to solve the equations of motion (|2.9|) - (|2.13|) 
by an analytical method, the boundary conditions (|3.1|) 
and (|3.2p are sufficient. However, Eqs. (|3.2p do not tell 
anything about the asymptotic behavior of the metric 
tensor at large p and thus a possibility of the gravity 
localization around the 3-branc remains unknown. The 
requirement of the gravity localization is equivalent to a 
finitcness of the four-dimensional Planck mass Mpi anc k, 
i.e., an inequality 



-^Planck 



dpM 2 {p)C{p) < oo 



(3.3) 



must be satisfied. Notice that solutions respecting the 
boundary conditions do not often satisfy the inequality 
(|3.3|) . Since the inequality requires a fine-tuning of pa- 
rameters for realizing the gravity localization, practically 
it works as the fifth boundary condition of the model. As 
we shall see below, imposing it on the model and con- 
sidering about a empty space-time with only the bulk 
cosmological constant, we can find the information of ge- 
ometry far from the origin. 



C. Vacuum solution of warp factors 

If a gravitational source in certain brane model is con- 
structed by a local topological defect, the all source terms 
in the Einstein equation will vanish for the region of 
p — ► oo. Namely an asymptotic behavior of geometry 
at infinity, which can not be determined by the bound- 
ary conditions, obeys vacuum solutions of the sourceless 
Einstein equation with the bulk cosmological constant. 

Combining the components of Einstein equation 
(|2.1ip - (|2.13[) . one can easily find equations for vacuum 
solutions, given by 



, 5 2 a fi 

m H — m H 

2 4 



0, I 



a(3 3m 
"im 2~' 



If the cosmological constant j3 is negative, then one ana- 
lytically obtains the solutions 

M{p) = M a e- Cp \l + ee 5cp \ 2 ^ 5 , 

C(p) = C e- Cp \ee 5cp - 1| • |1 + ee 5cp |- 3/5 , (3-4) 

where Mo, Co, e are arbitrary integral constants and the 
coefficient 



yj -a/3 /10 > 



(3.5) 



is a function of the model parameters. On the other hand, 
the case for the positive (3 was investigated in Ref . f55| . 
Therefore we shall concentrate the case of the negative 
bulk cosmological constant (3 at the present paper. Both 
cases are acceptable in seven dimensions [281 ] . 

Furthermore, in order to study the singularity struc- 
ture of the geometry, we investigate all curvature invari- 
ants for the metric tensor (|2.7|) . Explicit forms of the 
curvature invariants are 



K 



R 

kg 



A/TV 



M'C M' 2 ? M" 
ML " M 2 C M' 



RmnR 
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(kg? 



M' 2 C' 2 ni M' 3 C n M' 4 „M'CC" 
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TZabcdTZ 
,M' 2 C' 2 



ABCD 



RabcdR 



ABCD 



{kgf 



M' 4 M" 2 C" 2 



CabcdC 



ABCD 



CabcdC abcd 
(kg) 2 



12 [M' 2 C + M 2 C" - M {M'C + M"C)] ' 
5M 4 C 2 



(3.6) 



where 1Z,1Zmn,7Zabcd,Cabcd are defined as dimcn- 
sionlcss curvatures. The curvatures were already defined, 
e.g., in Ref. (36| . 

Inserting Eqs. (|3.4| ) into Eqs. (|3.6[ ). we study the geom- 
etry at infinity. The scalar and the Ricci curvature in- 
variants 

K = -30c 2 , K MN K MN = 150c 4 (3.7) 

are simply constants for any e, whereas the Riemann and 
the Weyl curvature invariants 



TIabcdTI 



P pABCD 

L-abcdL- 



ABCD 



60c 4 F(e), 
3840c 4 e 2 e 10c '' 



(1 + ee 5c P) 
are functions of e, where 

1 + iee 5c f + 70e 2 e 10cp + Ae 3 e 15cp + e 4 e 20cp 



(3.8) 



F(e) := 



(1 + ee 5c P) 4 



Thus we find that the scalar and the Ricci curvature in- 
variants are always constant and the regularity of the 
Riemann and the Weyl curvature invariants depend on e, 
for the above vacuum solutions. 

Next, in order to study whether the gravity can be 
localized and the geometry can be regularized on those 
solutions, we must analyze about the integral constant e. 
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For e = 0, the solutions become 
m (p) = t(p) = - 



and thus the warp factors 
M( P ) =M Q e- cp , 



(3.9) 



C(p) = C e~ 



(3.10) 

are exponentially decreasing. In the present case the Ric- 
mann and the Weyl curvature invariants are also con- 
stants similar to the others. Also each of the curvature 
invariants is equivalent to the one calculated from Ric- 
mann tensor for AdSjv space in Ref. [f56| . For AdS at space, 
each of the curvature tensors is represented in terms of 
the Gaussian curvature K, which are given by 

R = N(N-1)K, Rab = {N -l)KQ AB , 
Rabcd = K (GacGbd - QadQbc) ■ 



Since one can easily find K 



from TZ in Eas.([3~T|) 



in our six-dimensional model (namely N = 6), the other 
curvature invariants for AdSe can be also obtained sim- 
ply. Therefore the choice of e is desired one for the finitc- 
ness of the four-dimensional Planck mass (|3.3p ; since the 
asymptotic solutions can lead to a smooth AdSg geome- 
try far from the vortex string core and to localize gravity 
around the vortex string 2 . 

The integral constants e,A / foi£o have not been fixed 
for the above discussion. However, if a string-like defect 
is placed at the origin p = 0, the constants are no longer 
arbitrary and become functions of the model parameters 
a, /3, 7, namely e = e(a, (3, 7), and so on. The regular ge- 
ometry is achieved together with the gravity localization 
if the parameters lie on the surface e(a, /3, 7) = 0. There- 
fore we shall find numerical solutions with e(a, (3, 7) = 0. 



IV. ASYMPTOTIC SOLUTIONS 



In order to solve Eqs. (|2.9|) - (|2.13[) numerically, asymp- 
totic behaviors of the warp factors, scalar, and gauge 
fields in the vicinity of the origin as well as at large dis- 
tance of the core are mandatory [5?], [H, [|59[ . They are 
obtained by expanding the functions around the origin 
and approximating the equations at infinity. In this sec- 
tion, we also study about relations for a string tension 
since they may give useful informations to find the proper 
solutions. 



2 Other cases of e are as follows (see also Ref. |26|0 . If e > or 
e ^ — 1, we obtain exponentially growing solutions. For this case 
Eq.l l3.3l l diverges and therefore gravity can not be localized. (The 
case e = — 1 is somewhat specific since a singularity of Kasner 
type is developed in the origin.) If —1 < e < 0, the geometry has 
a singular point po and p < po should be required. In spite of the 
fact, the integral I l3.3t is finite. Thus this case has a possibility 
of localizing gravity if the singularity at po is resolved. Finally, 
let us note that if the bulk cosmological constant Af, is zero, 
solutions have a power-law behavior belonging to Kasner class. 
These solutions leave open only two possibilities but can not lead 
to localization of gravity. 



A. At the origin 

To examine behaviors at the origin, we start by ex- 
panding the warp factors together with the scalar and 
gauge fields as power series in p. Here we consider the 
case of n = 1. (For other n's one can estimate in a 
quite similar fashion.) Inserting the power series into the 
equations of motion (|2.9p - (|2.13p and requiring that the 
expanded equations obey the boundary conditions (|3.ip 
and (|3.2|) for a limit p — > 0, then one can get the asymp- 
totic solutions 

f( P ) ~TT + Ap, 

a(p) ~1 + Bp 2 , 

M( P )cl + -(-f3-2 1 + 2B 2 ) p\ 
C(p)cp-^[-P-2 1 + A 2 + 10B 2 ]p 3 , (4.1) 

where the coefficients A and B are two arbitrary con- 
stants. They can not clearly be determined by only lo- 
cally analyzing the equations of motion. We thus need in- 
formations of the results of numerical integration. Prac- 
tically, the constants are used to realize the boundary 
conditions at infinity. 

Furthermore, inserting the asymptotic solutions (|4.ip 
into Eqs. (|2.15p and (|3.6p . we find the asymptotic behav- 
iors of the energy-momentum tensor and the curvature 
invariants in the vicinity of the origin. We obtain the 
asymptotic forms around the origin 



r (p)^-2 1 ~^--2B 2 + O(p), 

T r (p)~-2 1 + 2B 2 + 0(p), 
Te( P )~~2 1 + 2B 2 + 0{p) 



(4.2) 



for the components of energy-momentum tensor and 
K ~ a [3 (/? + 2 7 ) + A 2 + 2B 2 ] + O(p), 



n AB n AB c ^ 



3 08 + 2 7 ) 2 + 2 (/? + 2 7 ) (A 2 + 2B 2 ) 



A 1 + \2A 2 B 2 + UB 



O(p), 



'Rabcd'R 



abcd 



O(p), 



CabcdC abcd ~ °-a 2 [A 2 



a 2 2(f3 + 2j-2B 2 ) 
+ ((3 + 2-1- A 2 - IQB 2 ) 

12B 2 -2{I3 + 2 1 )Y 
+ 0(p) (4.3) 
for the curvature invariants, respectively. 

B. At infinity 

Next we will consider asymptotic solutions for the two 
matter fields at large p. Let us assume that geometry is 



5 



regular at infinity as the case e = 0. Then asymptotic 
solutions of the warp factors are determined as Eqs. (|3.9p 

3 

From the boundary conditions for the matter fields 
(|3.I[) . asymptotic forms of them can be expressed by 

f( P ) = 6f(p), a(p) = Sa(p). 

Inserting the asymptotic forms and Eq. (|3.9p into the 
equation of motion for the gauge field (|2.10| . one eas- 
ily obtains 

(Sa)" - 3c{6a)' - 2(Sa) = 0. 

For p > 1 the solution a(p) can be described approxi- 
mately by 

8a(p)~e-v q = ~ h±yjl+^\. (4.4) 

Similarly, the linearized CP 1 field equation yields 

(Sf)" - 5c(Sf)' - 2 7 (5f) = 
which leads to the solution for the CP 1 field 

Sf(p)~e-», P=-*i(l±f^£)- («) 

In the asymptotic solutions (|4.4|) and (|4.5p . a minus sign 
should be chosen in order to obtain exponentially de- 
creasing behaviors. Since the asymptotic solutions must 
also satisfy the boundary conditions (|3 . 1 [) at infinity, the 
other exponentially growing solutions can never be al- 
lowed. 

For the completeness of our analysis, let us study the 
accuracy of our approximations more in detail. The equa- 
tion of motion for the CP 1 field (|2.9[) and the components 
of the energy-momentum tensor (|2.15p contain the terms 
of order C(l/£ 2 ), which exponentially diverge at large 
p. The linearization of Eq. (|2.9|) is justified if and only if 
q > 2c, i.e., 2 > a\(3\ which can be seen by easy calcula- 
tions. Also, as we shall see in Sec lIII C\ in order to drop 
the components of the energy-momentum tensor rather 
than the warp factors, we find that p > c and q > c, i.e., 
IO7 > 3a\(3\ and 5 > a\/3\ are required. Thus, especially 
for strong gravity limit, we need much thorough analyses 
for finding asymptotic solutions. 

C. Relationship for the string tension 

In four-dimensional cosmology, vortex-like topological 
defects often called cosmic string have been extensively 



studied. In particular, components of the string tension 
defined as 

/•OO 

Pm ■= / dpM i {p)C{p)T M {p) 
Jo 

contain many useful informations for the geometry [60l|. 
In the study of Abelian vortex in six dimensions [26| , a 
constraint is easily found by a reduction which we shall 
demonstrate in the below. Furthermore, the constraint 
plays an essential role in the stability analysis of vector 
mode fluctuations [Zi|[5(|. In our model a similar relation 
can be found but the usefulness is rather problematic: be- 
cause it contains a functional of f(p),a{p),M.(p),L(p). 
In the case of 'tHooft-Polyakov monopole in seven di- 
mensions [Hj], the situation is quite similar to ours. 

Consider two specific linear combinations of Einstein 
equation (|2.1ip - (|2.13p . namely 

1,1 r, a i \ a & 

m + 4m" + mi = —{t 8 + T r ) — , 

£' + e 2 + Ami = |(4r + r r - 3r e ) - ^. (4.6) 

Integrating Eqs. (|4.6p from zero to infinity, one easily find 
the relations 

lim M 3 {p c )M p (p c )£(p c ) 

= j(pe+Pr)-^-J° C M 4 £dp, (4.7) 
lim M i {p c )Cp{p c ) - I 

= ^(4p o +p r -3p )-^ J^^Cdp, (4.8) 

where the both right-hand-sides of them have already 
been taken to the limit p c — » 00. Let us note that we do 
not impose the localizing condition of gravity (|3.3p on the 
above formulations yet. In the limit p c — > 00, Eq. (|4.7p is 
the six-dimensional analogue of the relation determining 
the Tolman mass whereas Eq. (|4.8p is the generalization 
of the relation giving the deficit angle. Imposing the 
condition (|3.3p on Eqs. (|4.7p and (|4. 8|) and subtracting 

Eq. (|4.7p from Eq. (|4.8p . one obtains the following relation 

4 

Ho-pe = -— (4.9) 
a 

which must be hold where p r still remains undermined. 
In fact we do not use p r now. The condition (|4.9p is, 
however, only necessary but not sufficient in order to have 
solutions leading to localized gravity. 



3 The same geometry can be also realized in the case of t > or 
e < — 1, however both cases do clearly not realize the inequality 
(|3 - 3 1) for gravity localization. Thus we exclude these possibilities 
in subsequent consideration. 



4 If e 7^ 0, one can use other asymptotic solutions of the metric 
obtained by previous discussion for computation of the left-hand- 
sides of Eqs. H4.7ll an d l|4.8ll . Then this calculation gives a more 
general relation 12611 . 
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In order to get some informations of solutions at the 
origin, we directly integrate the left-hand-side of Eq. (|4.9j) 
and obtain 



Mo - Me 



A// 4 

dp^- {2v + d' 2 ) 



(4.10) 



Using the equation of motion for the gauge field (|2.10l) . 
the following relation is obtained 



M 4 d' 
£ 



M 4 



M 4 



f 



-^-v n = 2 — — I a — nsm' 
£ £ \ 2 



where Eqs. (|2.16[) are used at the last step above formu- 
lation. Inserting back it into Eq. (|4.10[) . the integral in 
Eq. (|4.10p can be rewritten as 



Mo - Me 




f 

a 2 + n (n - 2a) sin 2 ^ 



a' 2 



M 4 . 2 / , f M 4 a'a ^ ' 

I n(n — a) sin — \- 

£ v ' 2 \ £ 



2E 



f,a,M,£ 



M 4 a'a 



£ 



M 4 a'a 



(4.11) 



where for order we defined a functional of /,a, M,£ as 



f,a,M,£ 



M 



f 



dp—^-n (n — a) sin — 



For solutions the leading gravity localization, the bound- 
ary term at infinity in Eq. (|4.1ip exactly vanishes. More- 
over, using the boundary conditions at the origin and the 
fine-tuning relation (|4.9[) . we obtain 



a 



f,a,M,£ 



£ 



(4.12) 



According to the asymptotic solutions in the vicinity of 
the origin (|4.1[) . d(p) and £{p) are, respectively, a ~ 1 + 
Bp 2 and £ ~ 0{p) for the limit p — > 0. Finally, inserting 
these expansions into Eq. (|4.12[ ). we can find the relation 
with the model parameters a and B: 

B = -£- + s\f,a,M,c]. (4.13) 

As mentioned in the above discussion, the relation is not 
useful since it is given by the functional form. It may 
cause a serious problem especially on the stability anal- 
ysis describing later. 



V. NUMERICAL SOLUTIONS 
A. Method 

In this subsection the numerical strategy for our sys- 
tem will be outlined. An arbitrary second order dynam- 
ical system, in principle, can be expressed in terms of a 



first order system [57], [58|, [59J . Thus we transform our 
equations of motion (|2.9p - (|2.13p into first order differ- 
ential equations. By linearly combining the dynamical 
equations (|2.9[) - (|2.12[) . the following set of equations can 
be obtained 



y'i(p) 
y' 3 (p) 



y'e(p) 

y'r(p) 
y's(p) 



2/5 (p), 2/2 (M) = 2/6 (p), 

yz{p)yr{p)i y'±{p) = y&{p), 
2/8 (p)' 



4y 7 (p) H 
Vf{p) 



2/4 (m) 



2/5 (p) 



yi(p) 

42/7 (p) - 



2/9 (p) sin yi (p) 



2/8 (m) 
2/4 (p) 



2/6 (p) +V a (p), 



-\vi{p) + J [re( P )-(3], 
3 

-^y4(p)y 2 {p) - 3yi{p)y&{p) 

+ ^ V a(p) [4r (p)-3re(M)-/3] 



2/9 (P) = 0, 

in which we introduce functions yi{p) (i 
fined as 



1. 



(5.1) 
,9) de- 



2/1 (p) ■= f(p), 
Vi{p) := £{p), 
y 7 (p) := m(p), 



y 2 (p) := a(p), 
y 5 (p) := f'(p), 
y 8 (p) := £'(p), 



y 3 (p) :=M{p), 
2/6 (p) ■= a'(p), 
|te(p):=7- (5-2) 



Let us note that yg is actually not a function of p. 
The reason why we temporarily regard 7 as a function 
of p is as follows: the requirement of regular geometry 
e(a,/?, 7) = implies that 7 is a function of a and /3, 
i.e., 7 = 7(a,/3), which means that 7 must be uniquely 
determined once we shall give a and f3. In the sense, it 
is better to treat 7 as an eigenvalue of our model rather 
than a free parameter. We find that 7 becomes a constant 
once proper numerical inte gra tions shall be attained. 

The authors of Refs.p6l. l28l. HH have already pointed 
out that a 3-brane model described by certain field- 
theoretical Lagrangian almost contains enormous numer- 
ical difficulties. In language of a numerical analysis, the 
components of Einstein equation (|2.1ip - (|2.13p belong to 
the numerical class called the stiff equations which have 
exponentially growing solutions together with converging 
ones. Besides solving two-point boundary-value problem 
is another difficulty. 

The simple shooting method (SSM) is a very famous, 
relatively easy but an efficient method. Essentially SSM 
is a hybrid system of two solvers: one is for solving dif- 
ferential equations as the initial-value problem and the 
other is for matching boundary conditions. SSM is a 
simple and useful tool, however, necessarily must solve 
the initial-value problem in the interior, some instability 
may occur if the solutions strongly depend on the initial 
conditions. This is exactly the case that we must treat. 
Solving the equations of motion (|2.9p ~ (|2.13p by means of 
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SSM is thus a quite difficult task. As in the Refs. [2fl |28{ . 
we could also achieve a localizing gravity solution using 
SSM combined with the Down-hill simplex method 
At the present paper, we employ the parallel or multi- 
ple shooting method (PSM) [28|, [58|, [59| as the solver of 
the two-point boundary-value problem since PSM suc- 
cessfully dissolves the instability of SSM very well. 

The word of instability in our problem means that the 
functions released from a boundary-point tend to be di- 
vergent before arriving at the other. In PSM, in order 
to evade the instability of SSM, the integral region is di- 
vided into some subintervals where the running solutions 
is regular. After that, we apply SSM to each of them and 
continue the calculation until the functions converge at 
all junction and both boundary-points. Since PSM re- 
quires many initial conditions for all junction, called an 
initial trajectory, if one prepares an initial trajectory ap- 
proximate to proper solutions, it significantly saves the 
computing time. Thus the procedure to obtain solutions 
is as follows: (i) solve the equations by SSM with the 
Down-hill simplex method, (ii) use the solutions as an 
initial trajectory for PSM, (iii) solve by PSM to obtain 
new solutions, and (iv) repeat (ii)-(iii) until the conver- 
gence is attained. 

Another problem that we should care about is that 
the equations are so called the over-determined series. 
Our system is naturally constituted by the five equations 
(|2.9p - (|2.13[) but a number of the dynamical variables is 
only four. That means Eq. (|2.13[) works as a constraint 
of the system. Eq. (|2.13p is usually used to check con- 
vergence of the numerical integration. Also it may be 
a guide for finding the asymptotic solutions of the warp 
factors at infinity (in Sec lIII C| . However its information 
is never used in Eqs. (|5.ip . If one properly takes some lin- 
ear combination of the equations, all the five equations 
(|2.9p - (|2.13j) are involved in the numerical system. After 
that, 

4m ' - y mi = | (3r e - 5r r + 2/3) , 

g + I 2 + Ami = j (4t + r r - 'iT 6 - 2/3) 

or in numerical expression using Eq. (|5.2p 



' 1 1 1 1 1 


1 1 1 1 1 1 1 1 

— / 
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FIG. 1: A typical example of the solution leading grav- 
ity localization with the rescaled bulk radius p = \fkgr. 
The condition of shooting-parameters is as follow: n = 1, 
a = 0.9000000, [3 = -0.2000000, 7 = 0.488059791109, 
A = -1.60653214371 and B = -0.335537911342. 
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FIG. 2: An example of the energy-momentum tensor for the 
solution given by Fig[T] 



y'r(p) = f y7( P )^r\ + £ [3t,00 - Mp) + 2/3] , 
3 y 4 (p) 12 

2/8 0°) = - 4 2/70)2/sO) 

+ jVi(p) [4r (p) + Tr{fi) ~ 3r e (p) - 2/3] . (5.3) 

are found to significantly improve the numerical conver- 
gence because the information of fifth equation is taken 
into account. Thus the system can properly realize the 
asymptotic behavior at infinity. Henceforth we shall use 
Eqs. (|5.3p instead of the correspondences in Eq. (|5.ip . 



B. Results 

In this subsection, we present some typical examples 
in our numerical results. 

FigH] shows a typical numerical solution realizing the 
gravity localization around the 3-brane. Fig[2] and Fig[3] 
are the results of the energy-momentum tensor and the 
curvature invariants corresponding to FigfT] The results 
in Figf2] suggest that our solution describes a local topo- 
logical defect because all the component of the energy- 
momentum tensor vanish at large p. Furthermore, Fig[3] 
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2 4 6 8 10 12 14 

P 



FIG. 3: An example of the curvature invariants for the solu- 
tion given by Fig[T] From the considerations in Sec lHI Cl we 
interpret that the geometry achieves AdS6 where all curvature 
invariants are constant. 

certainly means that the geometry turns to be AdSg at 
large p. As already mentioned in Ref . [6l[ , the gravity lo- 
calization around a 3-brane, e.g., our configuration, may 
be a locally trivial incident of various cosmological events 
in the higher dimensional world. 

Also, Fig[4] shows the fine-tuning surface in the pa- 
rameter space (a, /?, 7) corresponding to solutions for the 
gravity localization. Contrary to the seven-dimensional 
case |28|, only /3 < 0, i.e., the negative bulk cosmological 
constant is allowed for localizing gravity solutions. The 
surface is growing for a — > and (3 — > — 00. In numer- 
ous studies [26|, HH, HH , the authors have insisted that the 
property of the solutions are dominated by the parameter 
c, i.e., Eq. (|3.5j) which relates to the asymptotic behavior 
of the warp factors at infinity and of course has much 
importance. Fig[4] tells us, however, that the dependence 
of the surface on the parameters a and j3 is complicated, 
thus more thorough analysis for physical implication of 
the parameter-dependence is required (see also FigJ5]and 
FigEJ) . In fact, 7 = j(a,/3) is the dimcnsionlcss coeffi- 
cient of the potential term and the variation of 7 brings 
remarkable change of the property of the solutions. 



VI. LINEAR STABILITY ANALYSIS 

In order to clarify the physical implication of our brane 
solutions, we proceed to linear stability analysis concern- 
ing the solutions. The stability of the solutions remains 
an open problem though it may be generally guaran- 
teed by its topology. In most of the previous works 
[47l . 0, Hi], HH , the source of gravity is well known 
Abelian vortex which is described by a scalar singlet and 
U(l) gauge field. The stability for the models in those 



FIG. 4: Parameter space for typical solutions which exactly 
localize gravity around the Maxwell gauged CP 1 brane with 
the winding number n = 1. Here a = \k is the dimensionless 
gravitational coupling constant, f3 = Ab/k 2 g is the dimension- 
less bulk cosmological constant and 7 = n/k 2 g is the dimen- 
sionless coupling constant of the potential like a mass term of 
a soliton. 

works is quasi-stable, which means that brane models 
described by topological solitons are not always engaged 
its stability. On the other hand, CP 1 model is written 
in terms of a scalar doublet (or a scalar triplet in terms 
of 0(3) sigma model) and is not based on gauge theory. 
Study of linear stability analysis for the present model is 
absent. It is thus worthwhile to examine its linear sta- 
bility afresh. 



A. Gauge-invariant framework 

The perturbed linearized Einstein equation is obtained 
by performing the first-order perturbation of the back- 
ground metric tensor. This procedure is straightforward, 
but one should care about the freedom of gauge, i.e., the 
choice of background coordinates before beginning the 
analysis. We adopt a gauge- invariant approach called 
the longitudinal gauge choice. The detailed discussion is 
given by Appendix [A] In this subsection we shall show 
only the outline of the analysis. 

The gauge-invariant fluctuations for the original metric 
tensor (|2.7|) can be written as 

(2M 2 U^ ME, MLT y \ 

SGmn = ME„ 2X LO , (6.1) 
\ MLFp LQ 2L 2 $ J 

where 

= V + r/^V (6.2) 

is the four-dimensional components. Also the gauge- 
invariant perturbations of the matter fields (|2.8[) is de- 
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FIG. 5: The behavior of the fields f{p),a(p),M{p) and £(p) for varying a at fixed f3 = -0.2000000. 



fined by 



Sz = (Exe- 

6A M = (e„ 



in I) 



d^0, O r , G e ) ■ 



(6.3) 



In the above equations, all the perturbed fields are func- 
tions of the full spatial coordinates x M = (x^ ,r,6). 

In order to derive the evolution equations of the fluc- 
tuations (|6.ip - (|6.3j) . it is convenient to introduce some 
perturbed quantities, e.g., the first order perturbation of 
Einstein tensor 



SG AB = SR AB - \ (GabSR + sg AB R) , 



where 



SR = G AB 6R A 



sg AB R AB 



(6.4) 
(6.5) 



is the perturbed Ricci scalar curvature. Furthermore, if 
one would like to find the explicit form of Eqs. (|6.4|) - (|6.5[) . 
one must calculate other perturbed quantities 

S gMN = -gMAgNB S g Ag ( g g) 

for the inverse metric tensor, 

STab = \<3 CD {^aSQbd + VbSGda - V D 5g AB ) 



(6.7) 



for the Christoffcl connection and 



SRicD = -doSTic + d c ST A D + ST E 



BD 0L EC 



BD EC 

8r% c ri D -T% c 6r£ D (6.8) 



for the Riemann tensor. In Eq. (|6.7p . Vjf denotes the co- 
variant derivative with respect to the metric tensor (|2.7[) . 
The evolution equations given by fluctuating the original 
equations of motion (|2.5p are thus obtained 



5G M n = X (STmn - Ai,5Gmn) , 
i 



6V M (T> M z a 



t a ?5z b = 0, 



8V M F MN = 

where 

STmn = 2kS 



S (A M V M z a ) 
N \z\ 2 ) 

+ iAkgS (z a d N z a * - z a *d N z a ) , (6.9) 



1 



V {M z a )*V N) z a 

+ (SQmn) Arane + GmN^C^ (6.10) 

is the perturbed energy-momentum tensor with the per- 



+ -6 [Q ab F M aFmb] 
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FIG. 6: The behavior of the fields f(p),a(p),M(p) and C{p) for varying /3 at fixed a = 1.000000. 



turbcd branc Lagrangian density 

<5£branc = ~kS [(%/)* V M Z°] - ^5 (F MN F MN ) 

-m[1-*(*W)]- (6.11) 

XimYn), in the first term of the right-hand-side of 
Eq. (|6.10[) . means as the symmetrization [38] defined by 



X(mYn) {XmYn + X N Y M ) . 

Since the evaluation of the evolution equations (|6.9p ~ 
(|6.11[) is straightforward but includes tedious algebra, we 
shall show the detailed descriptions in Appendix [Bj 

Note that the symmetrization in Eq. (|6.10[) naturally 
realizes the real gauge condition [5l|. The Si (i = 1,2) 
may be, in general, complex quantities, but we can ex- 
clude the imaginary part of Ei by applying certain U(l) 
gauge rotation. In mathematical point of view, the sym- 
metrization is equivalent to choosing the real part of a 
complex values. 



B. Tensor zero mode fluctuation 

For the tensor mode fluctuation, only the (fi, v) com- 
ponent of the perturbed Einstein equation is sufficient for 
the analysis. After a lengthy algebra, the final result of 
a evolution equation for the tensor mode fluctuation h^ v 
is given by 

^□v + K» + J* v + ( 4 ™ + K» = °> ( 6 - 12 ) 

where □ is the fiat four-dimensional d'Alembertian 

3 

n = v ^d^ = -d? + Y,df (6.13) 

i=l 

and 

dhiMir) s dlnLir) .„ 

rh := r -^-, I := -^- L . (6.14) 

dr dr 

In Eq. (|6.12p the prime and the over-dot denote the 
derivative with respect to the bulk radius r and the bulk 
angle 9, respectively. The first term in Eq. (|6.12|) can be 
n/iuy = m^h^L, for the massive graviton, since sat- 
isfies the four-dimensional Klein-Gordon equation. The 
term vanishes for the massless graviton. 
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Any physical quantity must exhibit axial symmetry 
since the geometry has same symmetry. We restrict 
to the massless (nih = 0) and the zero modes (p = 0), 
then 



K„{x M ) 



Eq. (|6.12p is fairly simplified to 



where 



is the dimensionless quantity and the prime, acting on 
variables with the tilde, denotes the derivative with re- 
spect to p. Also the functions m and £ in the equation 
arc defined by Eq. (|2.14j ). One easily sees that an arbi- 
trary constant is allowed as the solution of h^. Thus we 
employ h^ v (p) := h (constant). Since a normalizability 
of h^(p) is essentially equivalent to the integral (|3.3[) . 
it is automatically guaranteed if the background fields of 
perturbation satisfy the finiteness condition of the four- 
dimensional Planck mass (|3.3[) . 

Some authors in this field have arrived at similar con- 
clusions }49l . [50| . Since the tensor mode does not 
couple with the source term of gravity, properties of the 
localized massless graviton around a topological defect 
may be independent on the gravitational source in mod- 
els. 



C. Vector zero mode fluctuations 

Evolution equations for the vector mode fluctuations 
£ tl ^T il and <9 M are derived from the (p,, v), (p,r) and 
(p, 6) components of the Einstein equation and also the 
(/x) component of the equation of the U{1) gauge field. 
Here we introduce two new variables for convenience 



£p, ($J~ n) , ■ 



L 

M' 



(6.15) 



Similar replacements already have been used in Refs.|49l. 
[50j . From these variables, we obtain four equations 



4m + + = 0, 



1 



— D£ u + —€ u + ^ 



X 2a' 



M 2 ^ &L " qML 2 " 



On = 0, 



M 2 
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2 X fa', 



ML y g ^ + kv a O,)=0, 



(6.16) 
(6.17) 

(6.18) 



3m + # j 0'^ _ = 8kg0 ^ ( 6 .i9) 



We again concentrate on the lowest angular momen- 
tum eigenstatc upon the three vector modes. Then 
Eq. (|6~T7l) turns to be 



ip6 



p=0 



which means that D£ fl (x u ,r) = 0, therefore, for all inte- 
ger p, 



^□£ M (xV)e lpe =0. 



Then we conclude the vector mode £ ^ should be always 
the massless graviphoton. Inserting the condition into 
Eq. (|6.17|) . we obtain the important relation 



g M 2 



(6.20) 



which is valid for general angular momentum eigenstates 
since the vector mode £^ is the massless graviphoton. 
Also using the relation into Eq. (|6.18p . we obtain 



2x 1 
g ML 



4m — £) a' — kgv a 



O, 



0. 



For the last line we used the equation for the U(l) 
gauge field (|2.10[) . Consequently, the vector mode Tp. 
should be also the massless graviphoton for arbitrary an- 
gular momentum eigenstates. Those results mean, from 
Eqs. (|6.15|) and (|6.20[) . the vector mode M should be al- 
ways the massless gauge degree of freedom too. 

Next we insert Eq. (|6.20j) into Eq. (|6.19[) and define a 
dimensionless function with the rescalcd variable p 



kg ' 



we obtain the evolution equation for the rescaled vector 
mode 6> M (p) 

d"(p) + (2m + £)d'(p) 



a— -4]0»=O. (6.21) 



In the analysis of Abelian vortex [49j, |50j , all evolution 
equations for the vector mode fluctuation have been solv- 
able because the useful relation for the string tension was 
found. Unfortunately, in our case Eq. (|6.21[) is not ana- 
lytically solvable since the similar relation (|4. 9[) are less 
tractable for the analysis, as stated in Sec lIVI In this pa- 
per, we shall examine asymptotic solutions of O at both 
boundaries. At infinity, since we can linearize the evolu- 
tion equation (|6.2ip . we are able to find the asymptotic 
equation 

G"(p)-3c0'(p)- 89 M {p)=0 
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and the solution 
6> M (p) oc exp 



~ -1± 




(6.22) 



Also at the vicinity of the origin, using the asymptotic 
solutions (|4.ip for the background of perturbation and 
imposing a regularity of Eq. (|6.21|) at the origin, we find 
the asymptotic solution at the other boundary 

0» ~ 0<?) [1 + 2 (1 - aB 2 ) p 2 } + O (p 3 ) , (6.23) 

where is an arbitrary constant. From them, we spec- 
ulate that 6> M is almost localizing around the defect. 

For the vector mode £ ^ , we easily obtain the analytical 
solution from the constraint (|6.16|) . After rescaling the 
two vector modes with respect to p, i.e., 

:= £»{r), Puip) := ^>(r), 
we immediately achieve the equation 



£'(p) + (3m + £)£,(p) = 



with the solution 



e (0) 



where is an integration constant. For the vector 
mode J 7 from Eq. (|6.20p . we obtain 



M 



F*tfi) = '-7- \ 4 0) + 2a I dp 



C 



M 1 



where is also an integration constant. This solution, 
unfortunately, is a functional of a', 0^, and M.. 

Finally we examine normalizability of the three vector 
mode fluctuations. For £ M , T^, the integrals become 



d P M 4 C 



dp 



JO 2 

M 2 C 



(6.24) 



and 



dpM i C 



M 6 



dp 



/<°> + 2a / dp 



M 2 



(6.25) 



Clearly, they are not normalizablc if the background so- 
lutions lead localizing gravity around the 3-brane. More- 
over, the integrand of Eq. (|6.24p is certainly divergent at 
the both boundaries and that of Eq. (|6.25[) is also diver- 
gent at the origin. Thus we confirm that the two modes 
are not renormalizable without any information about a 
explicit form of the vector mode 0^. On the other hand, 
though we could not assure that the normalization inte- 
gral for 0^ is exactly finite, we confirm the integrand is 
finite at both boundaries from Eqs. (|6.22[) - (|6.23[) . In that 
sense, 0^ is the potentially normalizable vector mode 
fluctuation. 



D. Scalar zero mode fluctuations 

Since all perturbed equations contain the scalar mode 
fluctuations, we must deal with the whole evolution equa- 
tions (|6.9p . (Their explicit forms are given in Appendix 
iBl ) In Ref . [521] . in order to treat such complicated cou- 
pled system, the authors divide the whole system into 
many distinct subsets by using some redefinitions of the 
fields. As Ref.[5l|, however, we only have two subsets 
from Eqs. (|6.9[) if we restrict our analysis to the only 
lowest angular momentum eigenstates. One of subsets, 
called sector I, contains the scalar mode fluctuations 
(0, r , f2), and another subset, called sector II, contains 



1. Sector I: O, 6> r , 17 

This subset is constructed by the six evolution equa- 
tions of Eqs. (|6.9p ; the (p, 9) and (r, 6) components of the 
perturbed Einstein equation are 

Q' + 2 (m + i) n = (0 r - 0') - 2k x ^-0 (6.26) 



and 



2AI 2 



9 L 



UQ + [m -£ + Arh z - 1 - 3m£ Q 



= ^f2-k X ^9 r , (6.27) 
g L l L 

the (p) and (r) components of the perturbed U(l) gauge 
equation are 

0" -0' r + (2m + l) (0' -0 r ) = 8kg0, 



□ (0 r - 0') - jv a = 8kg0 r , (6.28) 

and two imaginary parts of the perturbed CP 1 equation 
in Appendix |Bj Since the present sector contains two 
constraint equations, we can perform to stability analysis 
of the sector with only three equations (|6.26[) - (|6.28[) 
Here we define a three dimcnsionlcss variables 

0{p):=0{r), r (p):=^P, Q{p) := 0(r), 

and also two new variables 

eiO>):=0»-0r(p), 02{p):=0(p)- 

Using these variables, we get the two equations from 
Eqs. (pr26l) and (jOg]) . 

Q'{p) + 2{m + l) Q{p) = (a'0!(p) + v a 9 2 (p)) , 



0[ (p) + (2m + €) 0i (p) = 802 (p) , 



(6.29) 
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FIG. 7: The magnitude of the function \ lQvC — v%\ for various 
parameter choices. 



and the two constraints from Eqs. (|6.27[) and (|6.28p . 



(6.30) 



in which we assumed all the three degrees of freedom are 
massless and used the relation 

m! -£' + Am 2 - I 2 - 3m£, 



a(r e -r ) = a ( — + — 



which is obtained from the components of original Ein- 
stein equation (|2.11|) and (|2.12|) and also the relation 
(|4.10p . In order to satisfy Eqs. ()6.30|) . there are two pos- 
sibilities: 

2v Va hi \ n 

~ a = sc or Q{p) = - 

The former condition obviously can not be achieved by 
the solutions of the numerical integration, e.g., shown in 
FigEl therefore we choose fi(p) = 0. 

From Eqs. (|6.30[) we find that the useful relation 



Inserting Eq. (|6.31|) into Eq. (|6.29p . we finally attain the 
evolution equation 



9'i{p) + (2m + I) 0' 2 (p) - W 2 {p) = 0. 



(6.32) 



Similarly to the above two subsections, this equation is 
not analytically solvable since m and i are the functions 
of p, but we can observe the asymptotic behaviors at both 
boundaries. At infinity the asymptotic solution for 2 is 



<9 2 (p) oc exp 



(6.33) 



From Eqs. (|6.31| and (|6.33p . the asymptotic solution for 
6>i is also 




0i (p) ex 02 (p). 



(6.34) 



The asymptotic solutions for 0\ and 02 are essentially 
the same as Eq. (|6.22p . On the other hand, we find the 
asymptotic solution at the origin by expanding the equa- 
tions and the functions as the same method in Sec lIVI 
The results are 



Oi(p) 

02(P) 



M 2 0} p + O( P 2 ), 

4 0) (l + 2p 2 )+0(p 3 ) 



(6.35) 



where ^ 2 °'' is an arbitrary constant. Therefore we suggest 
that 0i and 02 are potentially normalizable since the 
normalizability integrals for them are regular at the both 
boundaries. 



2. Sector II: V, X ,<P, E u E 2 ,Oe 

Analysis of the sector II is much more complicated than 
the sector I, since in the sector II we should treat the 
six degrees of freedom and the eight evolution equations 
which we shall sec in the below. Note that we again 
restrict the analysis for the lowest angular momentum 
eigenstates upon all the quantities. Firstly for the five 
components of the perturbed Einstein equation belonging 
to the present sector: the (p ^ v) component is 



d^d v (2$ + x + <p) = o, 



(6.36) 



e l {p) = e' 2 { P ). 



(6.31) the (p = v) component is 
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W" + <P" + 3 (4m + l) <?' - ( 3m + ^ I (A" - <?>') 



2 (3m' + 1' + 6m 2 + <? 2 + 3mi] + X = [a' (X + $) + &' e ] — /x% ( 1 — 2£i cos - + 2£ 2 sin - 
V / g C \ 2 2 



&X 



/' 2 „ , 2 ^ , «« 



2a 2 



— X + p $ + + /' sin /227J - -p- cos f2£\ - f cos /2£ 2 



, 2(n-a) 2 



L 2 



sin /2Z2 



, (6-37) 



the (r, r) component is 

jj^D (3<f + #) + 4 (3m + I) if + 4m'<P' = (<*'* + ©e) - MX ' 



/ / 
A" (1 — cos /) — Zi cos — + £2 sin — 



&X 



2«„ ... .,. /„, 2a 2 f ,,,/„, 2(n-o)* / 



p (* - X) + -^ 2 e - f sin J -S[ - cos L Sl + f cos J -£' 2 - 



sin -1^2 



2 X A 6 A\ (6.38) 



the (0, 9) component is 



M 2 



— □ (3& + X) + 4 + m (5^' - A") + (2m' + 5m 2 ) (<? - X) 



j!L(a'X-G' e )-nx< 



$ (1 - cos/) - 27! cos | + ^2 sin L 



r 2 



2 v " " y L 2 
and the (/i, r) component is 



(* - #) - ^ + /' sin J -S[ + cos J -£ x _ /' C os 



I 2 



sin — So 



2 x A b <P, (6.39) 



X a 



el / v-i /•/ „ / 



3<f + <£' - (m - <P - (3m + £J A = + fcx If sin J -£i - f cos J -E 2 

Next for the perturbed matter field equations, the real parts of the perturbed CP 1 field equations are 
nSi + £'{ + (4m + tj (e[ + Xf sin Q + Xf" sin | 



(6.40) 



2a 2 T 2a 



/ sin £ + ±—X + + —Og cos 



and 



^□£2 + ^2 + (4m + I) (z' 2 - Xf cos - Xf 



L 2 L 2 



"cos^ + ( 2 ^ + ^)/'co S | 



/ 


/ a 2 /i 


2 


\L 2 k 



£i=0 (6.41) 



/' 2 „ 2 (n - a) 2 ^ 2 (n - a) 

- — A H —<P -( 

2 L 2 L 2 



. f 
sm-- 



(n- a) M 
L 2 k 



£2 = 0, (6.42) 



and the only (6) component of the perturbed U(l) gauge field equation is 
-^□6> e + 6>^' + (4m - tj 0' e - a' (4<f' - X' - &) + 2Xv a = 8kg 

I 



/ / 
0$ — 2a£i cos — + 2 (n — a) Z" 2 sin — 



(6.43) 



The properties of each evolution equations are as fol- straints for all the fluctuations belonging to the sector II, 
lows: Eq. (|6.36p is a constraint for the three fluctua- and Eqs. (|6.37[) . (|6.39[) . and (|6.4ip - (|6.43p are dynamical 
tions of the geometry; also Eqs. (|6.38j) and (|6.40|) are con- evolution equations. Owing to the constraint (|6.36[) for 
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the geometry fluctuations, we are able to evaluate about 
the five degrees of freedom, i.e., X, S\, S 2 , and Qg. 

If we consider the special situation where all five fluctu- 
ations are massless, we obtain the asymptotic equations 
at infinity from Eqs. (|Qg|) . and (|^IT|) - f05|) . 

2&(p) X'{p) 5cX(p) = ^ (l - 2Ei{p)) , 
§'(p) - X'(p) + AcX{p) = (6.44) 

for the fluctuations of the geometry and 

G' e \p) - 3cO' e (p) = 80 e (p), 
Z'{(p)-5cZ' 1 (p)+jZ 1 (p)=0, 
^(p)-5cE' 2 (p)- 1 S 2 (p)=0 (6.45) 

for the fluctuations of the matter fields, where we again 
use the variables rescaled with respect to p. Consider the 
fluctuation of the constraint for the CP 1 doublet \z\ 2 = 1, 
that is, 



Si cos — + £2 sin — 



S\z\ 2 = 5 (z* a z a ) = Sz* a z a + z* a 5z a 

5(1) = 0, 

we find the constraint for the perturbed CP 1 fields 
S 1 (p) = £ 2 (p)tan^. 

£\(p) goes to zero at infinity since it must hold the 
boundary conditions for the backgrounds (|3.ip . Hence 
we determine the asymptotic solutions for Eqs. (|6.45D 



O e {p) <x exp 
S 2 (p) oc exp 



— 1 ± l/l 



32 
9? 



— 1±\ 1 



87 
25c2 



(6.46) 



Also, we get the asymptotic solutions of the fluctuations 
for the geometry from the above results and Eqs. (|6.36p 
and (jnUll): 

*(p)^ + x (oo) e 13 ^, 

4c 

31 



S(p) 



^L + 2a 1 p~^x ioo) e^, (6.47) 



where ^(°°) ; and <fr(°°> are arbitrary integral con- 
stants which satisfy Eq . ([OB")) . Egs- fOB]) and ([6^47]) 
clearly show that the three scalar mode fluctuations of 
the geometry diverge at infinity even though the other 
three fluctuations of the matters are regular. 



On the other hand, we investigate asymptotic behav- 
iors at the origin as the same method in Sec lIVl Conse- 
quently we get the asymptotic solutions as follows: 



Hp) - - g (7 
X{p) ~ x ( V 



0(P 3 ), 



o( P 3 



e e ( P )^d e 0) p 2 + o(p 3 ), 



+ o( P 5 ), 



(6.48) 



where , ; an d are arbitrary constants. For the 
rest of degrees <P(p), we similarly evaluate by plugging the 
asymptotic solutions into the constraint (|6.36|) . 

Therefore, from Eqs. (|6.46[) - (|6.48[) . we conclude that 
the geometry parts of the fluctuations belonging to the 
sector II are not normalizable degrees whereas the mat- 
ter parts of them are potentially normalizable ones. Of 
course, we should be practically required some numeri- 
cal calculations in order to find exact solutions for the 
evolution equations. 



VII. CONCLUSION 

In this paper we have investigated a new brane model 
in six dimensions, constructed the thick brane solution 
by the Maxwell gauged CP 1 model. The origin of this 
model is naturally obtained by projecting 0(3) sigma 
model onto a complex space and exchanging a P(l) 
gauge term with the native composite connection. As 
we have shown, our model can realize localizing gravity 
around the 3-brane with certain parameter space of the 
model and attain a finite four-dimensional Planck mass 
fEq. (|3.3p ). Our results clearly suggest the possibility of a 
new variety of braneworlds based on classical field theory 
because there are many variants for 0(3) sigma model 

[u, m, m, m, m m m, m m m (see also Ref . m ) . 

Besides, relevant previous works are almost constructed 
from the classical solutions in gauge theory; in partic- 
ular Abelian-Higgs (AH) vortex in six dimensions [2(| . 
By using the different model from AH vortex, we have 
showed that a solitonic nature of vortices can attain grav- 
ity localization around the brane core in codimension- 
2 braneworld models. Since such mechanism using any 
topological soliton models does not always work in higher 
codimensional models (l9| , it is significant that our model 
practically obtain similar results to AH vortex . 

Our another aim at the present paper is to analyze 
the linear stability of the model, for all fluctuating fields 
around the background classical solution. We have con- 
centrated to study only their zero modes here. For the 
fluctuations of the geometry, the tensor zero mode is 
always localized around the origin, on the other hand 
neither the vector nor scalar modes are localized. The 
fluctuations of the matter fields are potentially localized 
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since they are regular at infinity and near the origin. 
Here, "Potentially" means that we could not find any an- 
alytical solutions. In order to determine whether these 
fluctuations of matters are true stable or not, we need 
vast numerical calculations for all the perturbed equa- 
tions of motion presented in Appendix [Bj We therefore 
conclude the brane described by Maxwell CP 1 model is 
a quasi-stable configuration. 

For the linear stability analysis on the zero mode fluc- 
tuations, other thick brane models in six dimensions 
[49l . 1501 15^ | have obtained similar results to us. These 
codimcnsion-2 models inclu ding ours have a serious draw- 
back as pointed out in Ref. [5l| . One of aims in our pa- 
per is to study such problem by using other model with 
different topology from AH vortex, but we finally en- 
counter a similar difficulty. However, we believe that 



thick brancs include rich physical implications than thin 
ones, thus we should tackle with the problem. In recent 
years, the many attempts to resolve the problem of the 
quasi-stability by introducing other degrees of freedom. 
We are also currently studying such quasi-stability prob- 
lem about our solutions and the results will be reported 
in near future. 
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APPENDIX A: LONGITUDINAL GAUGE SYSTEM 



Mathematically, the problem of describing evolution equations for small perturbations in general relativity is equiv- 
alent for solving the Einstein equation linearized about an expanding background metric tensor. This procedure 
is straightforward but there are complicated issues concerning the freedom of gauge, i.e., the choice of background 
coordinates. In this appendix we derive the gauge- invariant fluctuation for all fields (|6 . 1 1) - (|€3 . 3[) in Scc lVII bv imposing 
general fluctuations upon the longitudinal gauge conditions. 

First of all, we introduce the first-order general fluctuation of the background metric tensor (|2.7p given by 



sg 



MN 



S^Gmn + S^Gmn + S^G 



MN 



MV, 
V MLWu 



2* 
Llu 



Llo 
2L 2 d 



(Al) 



where 



Uu,v — h. 



with the divergenceless and the traceless constraints 



cU£ = 0, M = 0; d,P L = 0, 8^ = 0, dJC> i = 0. 



(A2) 



In above fluctuations, is a tensor mode, X M , J^, and /C M are vector modes, and finally -0, 0, V , Q, Xi 4>i an d w are 
scalar modes, which depend on all the spatial coordinates x M = (x^ 1 , r, 8). 

As we mentioned before, any gauge-invariant quantities are required to be invariant with respect to the choice of 
background coordinates. Here an infinitesimal coordinates transformations are defined by 



A/ 



+ t M (x N ) 



(A3) 



with 



/■Mr N\ I 71 r 2 11 r T 2 0\ 



where the local infinitesimal parameter e M depends on x M . Also we take the four-dimensional sector of the transfor- 
mations as follows: 

e M = ^ + C- (A4) 
If the perturbed metric tensor (|A1|) is transformed SGmn under the transformations (|A3|1 . which can be written by 

SGmn = SGmn - £^Gmn- (A5) 
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where £j is the Lie derivative defined as 



M-. 



where Vm means the covariant derivative with respect to the metric tensor. When one calculate the Lie covariant 
derivative involved in the perturbed quantities and others, one may use the original metric tensor (|2.7| and the original 
Christoffel connections 



A// 7 T 1 



-LL'. 



Inserting the perturbed metric tensor (|A1[) into Eq. (|A5p . we obtain the explicit forms of SGmn- Taking some linear 
combination of them, we can obtain following forms of the scalar functions 



& := if) - M' [V - MO' 

m(v- MO' 

ML' 



X:=x 
3:=d>- 



L 



V - MO' - 



M 



M ~ 



— \Q--0 



Q := ui 

L 



M(V - MO' 



L 



M 



L 



M ~ 



— [Q-—0 



L 



(A6) 



where the prime and the over-dot denote the derivative with respect to the bulk radius r and the bulk angle 9, 
respectively. Similarly two gauge-invariant vector functions 



M ~ 

■= Jii - MI^, := /C p - —Ifj, 



and one gauge-invariant tensor function 



are obtained. 

On the other hand, for the matter fields, the original CP 1 and U(l) gauge fields (|2.8p are fluctuated by 



Sz = {SzuSztf = (a 1 (x M )e- m ^,a 2 (x M )) 1 ; 5 Am = (^(x iv ), tf r (x iv ), d e (x ly )) 



M\\T. 



(A7) 



(A8) 



(A9) 



with the decoupled vector component 



where 6> M and $ are a vector and a scalar fluctuation for the U(l) gauge field, respectively. We obtain transforms of 
the matter fields under Eq. (|A3j) in the same way as the geometry. Hence we can defined six gauge-invariant scalar 
functions given by 



Et := 5-1 + — M (P - MO'^j sin - 
E 2 ■= 5-2 - (P - MO') cos | 



M 

T 



M - 



L 



O cos 



(A10) 



for the CP 1 field and 



M 



0u = 0„, 0:=d-A e —lQ-—0\ 1 



r \= § r - A 6 

0s := $o - A 6 



L 



M - 









It 


(«- 


n 












It 




n 









L 



-A 



ML' 



L 2 

- A'nM ( V - MO' 



(All) 
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for the gauge field where the vector mode (9 M is automatically gauge- invariant variable like h^. 

Finally we take the longitudinal (or Conformal Newtonian) gauge conditions [371, where 0,V, Q,ifi are to be zero 
after the infinitesimal coordinates transformations (|A3[) . These gauges can be written in terms of 0,V, Q,^ as 

e = 0, e r = M(V-MO'), e e = ^ \Q - ^-O^j , 

and = 1^. Then the gauge-invariant fluctuations of the geometry (|A6[) - (|A8[) and of the matter fields (j A10|) and 
(jAlip are exactly equivalent to the original one under the transformations (|A3|) . In our perturbative calculation, 
therefore the independent gauge degrees of freedom are X, <l>, Q 7 T^, h^ Vl E±, £2, (9^, O, ri and Og. 
Also we omit the tilde for convenience in Sec lVIl 



APPENDIX B: PERTURBED QUANTITIES 



The explicit evaluation of the perturbed equations of motion (|6.9[) is proceeded here. We shall show the clear forms 
of the perturbed Christoffcl connection, Einstein tensor, energy-momentum tensor, and equations for the matter fields 
in the below. 



1. Christoffel connection 



Inserting the explicit form of the gauge-invariant perturbed metric tensor (|6.1[) into Eq. (|6.7p . one can get all 
components of the perturbed Christoffel connection after lengthly calculations. The (/i, M, TV) components are 



° re ~2M + 2M 



L 

2AP 



the (r, M, N) components are 

ST r al3 = —M 2 (U' a/3 + 2mU a(i ) + M8 {a S 0) + 2M 2 i la pmX, 
T' a + ( m + 1) T a 



_M- ML 
dL °e - 2 a 2 

5r; e = x-l'o, 



SI* r =U£ + — {d a £»-d»£ a ) l 



ST r ar = -M'Sa + d a X, 

sr r rr = x 1 , 

STffg = Lfi — L 2 



<P' + 2£(<P- X) 



(Bl) 



(B2) 



and the (8, M, N) components are 
M 2 



ST 



T 



ST, 



M 

2L 2 
1 



f 



M r 

2L 



T' 



2L 



d a f2, 



;X 



1 



L 2 " ' L 
STL = <P + L'Q. 



n 1 + in 



(B3) 



2. Einstein tensor 



Using the explicit forms of the perturbed Christoffel connections which we obtained, one can also evaluate the 
perturbed Riemann tensor (|6.8p and Ricci scalar (|6.5|) . Inserting these perturbed quantities into Eq. (|6.4p . one can 
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achieve six components of the perturbed Einstein tensor. The four-dimensional component (//, v) 

M 2 ■■ r / a / 

5G^ = -□ V - ttV " M ' [K» + ( 4 ™ + £ ) tip* ~ 2 V ( 3r ™' + £ ' + 6 ™ 2 +£2 + 3M 



d {tl £' v) + (3m + £J + — d^T v) - (d^d u - (2<F + Af + <P) 

M 2 r\ tiv | (3# + A?) + 3^" + <2>" + 3 (4m + f) - (3m + /) (A" - <2>') 

2 (3m' + ? + 6m 2 + I 2 + 3ml) (X - <P) - j 



Sl' + [Zrh + i ) h 



has all mode fluctuations of the perturbed quantities. The three extra-dimensional components (r, r), 



<5G rr = -^□(3<F + #) + - (y#-mJ7) +4(3m + l)l? / + 4m#', 



SGee = jpP ( W + x ) + L 2 [^" + 4 ™ (5#' - A") + 4 (2m' + 5m 2 ) (<P - X)} 



SG 



L 



2AP 



Uft + 2L (2m' + 5m 2 ) [2-4 



4>' + [rh-i] & -mX 



have the scalar mode fluctuations only. Finally, the rest of component (/i, r) and (/x, 
8G^ r 



L D£ u - —5-4 + M ( 3m' + I' + 6m 2 + I 2 + 3ml) £ u + A 



2M 



2L 2 



2L 



I X - —SI 



L 



M r . 



p 2M * 2 



£'^ + [bm-l)£ ll 



ML 



(B4) 

),and (r,9) 

(B5) 

(B6) 
(B7) 



(B8) 



I + Urn + tj - (jm! + £' + 17m 2 + 2l 2 + mt\ 



n' + 2 ( m + i ) n 



-w - x 



(B9) 



have the vector and scalar mode fluctuations. 



3. Energy-momentum tensor 



In order to evaluate the perturbed energy- momentum tensor (|6.I0[) . one need the explicit form of the perturbed 
matter Lagrangian density (|6.1I[) . Calculating all the fluctuation terms one can obtain 



SC 



branc 



f' 2 V a „ 2V 



J —X + -^Og + —$ + f sin J -S[ - 2jz cos J -Zy - /' cos J -E' 2 - 2 

a 1 (X + <P) + O'g - O r 



2" 1 ~L 2 ~~~2' 
gL 2 



L 2 S[U 2 S2 



f f 
[i [ 1 — 2Si cos — + 2U 2 sin — 



(BIO) 



Hence all components of the perturbed energy-momentum tensor can be given as follows. The four-dimensional 
diagonal component (/i, v) 



1 a! 



ST^ = —Mr]^ < -—y a '{V-X-$) + O r - 0' e 



f f 

1 + 2V (1 - cos /) - 2Ei cos J - + 2S 2 sin J - 



r il(*-X) + f 2 {*-<P)-^ 2 -O e -f sin ljJ[ + 2^ cos f -£ l + f cos f -E> 2 + sin f -S 2 



- 2M 2 h llv 



12 v \ 1 a' 2 



(Bll) 
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has the tensor and the scalar mode fluctuations, even though the corresponding component of the perturbed Einstein 
tensor contains the all mode. The two extra-dimensional diagonal components (r, r) and (9, 6) 



2v 



a'<P + 0' -0 r ) -n 



f f 

1 + 2X (1 - cos /) - 2£i cos - + 2E 2 sin - 



f r-i/ 



— (<P-X) + ^0 e - f sin J -S[ - cos J -S x + f cos J -Z 2 



f o ("a ~ a) . / 



L 2 



sin — S2 
2 



ST e , 



1 



9 



a' ( a'X + O r - 0' fi 



f f 

1 + 2<£ (1 - cos /) - 2Z 1 ! cos J - + 2£ 2 sin J - 



f' 2 



f 1^' 



- — (X — $) — -^Og + f sin Ls' + 2— cos - f cos J S 
2 v ; L2 t j 2 1 L 2 2 2 2 



/ o ("a - a) . / 



sin — 

L 2 2 2 



have the only scalar mode fluctuations. The two off-diagonal mixing components (/i, r) and (/x, 



g 1/ 
1 



- kdf, ( /' sin -Si- f cos -Z2 



<ff> = ML^C hxme + -a' [d„ (0 r - O') - ©;] - fcu (<9 M <9 + ©„) 

have the vector and the scalar mode fluctuations. Finally, the rest (r, 9) component 

la' 2 / / • f ■ 

ST r g = LJ?£br a nc + _ fc (/' Sm T^ 1 ~ ^ C ° S ? ^ 2 + ^ a<9r 

has the only scalar mode fluctuations. 



(B12) 



(B13) 

(B14) 
(B15) 

(B16) 



4. Matter field equations 



For the CP 1 and U(l) gauge fields, one can obtain the perturbed equations of motion by inserting the perturbed 
fields (JO]) into Eq. fOj) . 

For a = 1 component of the evolution equation for the CP 1 field, the real part is 



M 2 



as 1 + s'l + + Urh + t )[S[ + Xf 



,/-:_/ 



Xf" 



f r\ . f 
— il sm ■ 
2L 



29' 



#-X'\ , , f , (f 
f sm - -' 



til 



-X 



2a' 
T 2 



(p 



2a 



Og cos ■ 



a 

L 2 



27i = 0, (B17) 



and the imaginary part is 



('no 

\M 2 



1 



L 2 



f 



f 



U0 + 0' + 4*06 ) cos i- + (4m + I ) (0 r + ^n) cos J - - [AW + X - $) cos J - 



f 



aft' + a Q — aQl) cos 



L 2 
a 
L 



)f si nt-^A = 0. (B18) 



For a = 2 component, the real part is 



1 



M 2 



29' 



L 2 
& - X' 



UE 2 + E' 2 ' + -prSi + Urh + £) £ 2 — Xf cos 



f 



Xf" 



r r 



2L 



n cos : 



J 2 



f 2 
—X 
2 



2 (n a - a) _^ _ 2 (n a - a) Q 



L- 



L 2 



. f 
sm — 
2 



L 2 



S 2 = 0, (B19) 



and the imaginary part is 



(J- 



D0 + 0L 



1 

L2' 



Og sin 



/ 



a'n+(n a -a)(n£- 



L 

. / 
sm — 
2 



H sin 



/ 



Or - 



49 + X - <P 



n f cos 



n a - a . J 

7. — sm — 

L 2 2 

/ 2 (n a - a) 

2 L 2 



S 2 = 0. (B20) 
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There are the only scalar mode fluctuations in Eq. (|B17|) - (|B20[) . 

Similarly to the CP 1 field, one can obtain the evolution equations for the U(l) gauge field. The (fi) component of 
this equation has the vector mode sector 



1 



1 



M 2 M L 2 
and the scalar mode sector 



—:Ue u + + —0 U + (2m + i )9 ' - —a' I -£ u - 



L L 



T'„ + (i-m) T, 



(B21) 



0" -G' r + (2m + I) (0' - r ) + [0 - 0g ) = 8kg0. 



In contrast to the component (/x), the (r) component 



1 



1 



jj^U (0 r + [Or -0' e )+jj{W-X-<P)- -v a = 8kg0 



and the (0) component 



(B22) 



(B23) 



-^2 □ {S e -0)+ O'g' -0' r + (4m - tj (e' 6 - 6> r ) - a 1 (AW 1 - X' - $') + 2Xv a 



= 8kg 



00 - 2aS\ cos — + 2 (n a - a) S 2 sin — 



(B24) 



have the only scalar mode fluctuations. 
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